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ABSTRACT 

We obtain an exact, closed-form expression for the time-dependent Green's function 
solution to the Kompanccts equation. The result, which is expressed as the integral 
of a product of two Whittaker functions, describes the evolution in energy space of a 
photon distribution that is initially monoenergetic. Effects of spatial transport within 
a homogeneous scattering cloud are also included within the formalism. The Kom- 
paneets equation that we solve includes both the recoil and energy diffusion terms, 
and therefore our solution for the Green's function approaches the Wien spectrum at 
large times. This was not the case with earlier analytical solutions that neglected the 
recoil term and were therefore applicable only in the soft-photon limit. We show that 
the Green's function can be used to generate all of the previously known steady-state 
and time-dependent solutions to the Kompaneets equation. The new solution allows 
the direct determination of the spectrum, without the need to numerically solve the 
partial differential equation. It is therefore much more convenient for data analysis 
purposes. Based upon the Green's function, we derive a new, exact solution for the 
variation of the inverse-Compton temperature of an initially monoenergetic photon 
distribution. Furthermore, we also obtain a new time-dependent solution for the pho- 
ton distribution resulting from the reprocessing of an optically thin bremsstrahlung 
initial spectrum with a low-energy cutoff. Unlike the previously known solution for 
bremsstrahlung injection, the new solution possesses a finite photon number density, 
and therefore it displays proper equilibration to a Wien spectrum at large times. The 
relevance of our results for the interpretation of emission from variable X-ray sources 
is discussed, with particular attention to the production of hard X-ray time lags, and 
the Compton broadening of narrow features such as iron lines. 

Key words: radiation mechanisms: thermal — radiative transfer: line profiles - 
plasmas — galaxies: active — cosmology: early universe — methods: analytical — 
X-rays: general. 



1 INTRODUCTION 

In hot, radiation-dominated plasmas, the primary interaction between photons and electrons occurs via Compton scattering. 
This process consequently plays a central role in the formation of X-ray spectra in a variety of sources, including active galaxies, 
low-mass X-ray binaries, and the early universe. The repeated scattering of radiation off free electrons (Comptonization) can 
have a profound effect on both the X-ray spectrum and the Fourier structure of the observed emission in the time domain. 
The effect of multiple Compton scattering depends on the velocity distribution of the electrons; random and systematic 
motions of the electrons produce "thermal" and "bulk" Comptonization, respectively. Thermal Comptonization can explain 
the formation of the power-law tails observed in the X-ray spectra from active galaxies and low-mass X-ray binaries, as well as 
the development of time lags between the hard and soft photons (Shapiro, Lightman, & Eardley 1976; Sunyaev & Titarchuk 
1980; van der Klis et al. 1987; Stollman et al. 1987). This process also governs the evolution of the cosmic background radiation 
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before recombination (Sunyaev & Zeldovich 1970; Zeldovich & Sunyaev 1969; Illarionov & Sunyaev 1975a, 1975b), and it is 
a necessary ingredient in models for the production of the modified blackbody continuum emission observed during X-ray 
bursts (Sunyaev & Titarchuk 1980; Titarchuk 1988). Thermal Comptonization also seems to contribute to the broadening of 
the Ka iron lines observed in a number of sources (Reynolds & Wilms 2000; Wang, Zhou, & Wang 1999; Ross, Fabian, & 
Young 1999). In addition to the thermal process, bulk Comptonization can also influence the spectral formation when the 
electrons are converging on average, as for example in a radiation-dominated shock or a quasi-spherical accretion flow (Payne 
& Blandford 1981; Lyubarskii & Sunyaev 1982; Becker 1988; Colpi 1988; Titarchuk, Mastichiadis, & Kylafis 1997; Titarchuk 
& Zannias 1998; Laurent & Titarchuk 2001; Titarchuk & Shrader 2002). 

The effect of thermal Comptonization on the photon distribution is described by the well-known Kompaneets partial 
differential equation, published in 1957. This equation, based on a Fokker-Planck formalism, describes the effect of multiple 
Compton scattering on the photon distribution when the electrons are moving nonrelativistically and the average fractional 
energy change per scattering is small. The Kompaneets equation has provided the foundation for most of the theoretical studies 
of thermal Comptonization in astrophysical environments, which have focused primarily on the steady-state reprocessing of 
monoenergectic radiation injected continuously into a nonrelativistic plasma with time-independent properties (e.g., Katz 
1976; Sunyaev & Titarchuk 1980). These studies have provided useful insight into the formation of the ubiquitous power-law 
spectra observed in the quiescent emission from active galactic nuclei (AGNs) and galactic black-hole candidates. Titarchuk 
& Lyubarskij (1995) generalized these results by solved the stationary form of the Boltzmann kinetic equation to demonstrate 
that power-law spectra can also result from Comptonization in a cloud of relativistic electrons. Despite the success of the 
steady-state models in helping us to understand the formation of the quiescent X-ray spectra, they cannot be used to follow 
the spectral and temporal evolution of rapid transients. This is a concern because rapid variability is a characteristic feature of 
many classes of X-ray sources. Indeed, the rapid variability itself contains detailed information about the geometry and spatial 
structure of the inner region of the accretion flow. Analysis of the "Compton reverberations" associated with the impulsive 
injection of soft photons can provide perhaps the most direct probe of the inner region (Reynolds et al. 1999; Ulrich 2000). 
Specific examples of rapid variability include the strong X-ray flares observed from the Seyfert galaxies NGC 4151 (Tananbaum 
et al. 1978; Lawrence 1980) and NGC 6814 (Konig et al. 1997; Mittaz & Branduardi-Raymont 1989) in which the intensity 
increases by a factor of 5-10 in < 1000 s. Many other active galaxies also display significant variability (see, e.g., Iwasawa et 
al. 2000; Merloni & Fabian 2001), and flickering and bimodal spectral behavior on timescales of seconds involving sustantial 
fractions of the X-ray flux have been observed from Cyg X-l (Tanaka 1989; Zdziarski et al. 2002), GX 339-4 (Kitamoto 1989; 
Miyamoto et al. 1991), and other galactic black-hole candidates (Mereghetti 1993; Smith, Heindl, & Swank 2002). 

Observations of rapid variability during strong X-ray transients in AGNs and galactic black-hole candidates suggest that 
a time-dependent model for the spectral formation may be required in order to understand the physical processes occurring 
in the source plasma. Surprisingly, the time-dependent Green's function associated with the Kompaneets equation has never 
been obtained in closed form. Time-dependent analytical solutions previously discussed in the literature have focused on 
certain special cases, or have adopted simplifying approximations. For example, Chapline and Stevens (1973) examined the 
Comptonization of a bremsstrahlung initial spectrum in a plasma with steady properties. However, the time-dependent 
analytical solution they obtain is unphysical in the sense that full equilibration to a Wien spectrum never occurs because the 
initial spectrum contains an infinite number of low-energy photons. Zeldovich & Sunyaev (1969) and Payne (1980) derived the 
Green's function for time-dependent thermal Comptonization under the assumption that e <C kT e , where T e is the electron 
temperature and e is the photon energy. Under this restriction, the effect of electron recoil can be ignored, and the transport 
equation admits a relatively simple solution for the Green's function. However, this "soft-photon" Green's function is not 
valid when photons remain in the plasma long enough to upscatter to energies comparable to the electron thermal energy, 
which cannot be ruled out in the observed transients. Futher discussion of this point is provided in § 6.1. 

More general results have been obtained using numerical simulations based on a modified form of the Kompaneets equation 
that includes an additional term describing spatial diffusion. For example, Bottcher & Liang (1998) simulated numerically the 
Compton upscattering of flares of soft radiation in clouds with various geometries, and Malzac & Jourdain (2000) utilized a 
Monte-Carlo code to model time-dependent Comptonizing flares in a corona with a self-consistently determined temperature. 
The utilization of partial differential equation solvers or Monte-Carlo simulations to analyze the transport equation is a 
complex procedure that is inconvenient from the point of view of X-ray data analysis. Furthermore, the results obtained for 
the spectra often depend on boundary conditions in the energy coordinate that are not known with precision a priori. In order 
to limit the effect of the imprecisely known boundary conditions on the computed spectrum, one is usually forced to extend 
the computational domain well beyond the region of interest. The associated decrease in computational efficiency reduces 
the utility of simulations based on numerical integration of the transport equation. Moreover, even when accurate numerical 
solutions can be obtained, it is often difficult to extract simple analytical estimates from them due to the complexity of the 
simulations. 

Our increasing capability to make observations with high temporal and spectral resolution presents us with an inter- 
esting theoretical challenge, which is to obtain a closed-form expression for the X-ray spectrum resulting from the thermal 
Comptonization of initially monoenergetic photons in an ionized plasma. The availability of such a solution would be of great 
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importance for our understanding of spectral formation during rapid X-ray transients, and particularly when the Compton 
broadening of narrow features such as iron lines is of interest. As a first step, in this paper we analyze the transport equation 
for thermal Comptonization and derive the time-dependent Green's function for the case of a homogeneous plasma with 
steady properties. Our goal is to obtain the "complete" Green's function, including the effect of electron recoil as well as 
the diffusion of photons in the energy space. The formalism also incorporates a treatment of spatial transport inside the 
cloud. The resulting solution can be used to obtain a comprehensive understanding of the spectral evolution due to thermal 
Comptonization in a plasma with any shape and with any temporal and spatial distribution of embedded photon sources. 
We believe that the insights gained using this new solution will be valuable for the interpretation of current and future X-ray 
data from active galaxies and low-mass X-ray binaries. 

The remainder of the paper is organized as follows. In § 2 we introduce and discuss the fundamental transport equation 
governing the propagation of photons in X-ray emitting plasma clouds. In § 3 the transport equation is solved to obtain the 
formal solution for the time-dependent Green's function in the energy domain as an integral of a product of two Whittaker 
functions. The properties of the Green's function are explored in § 4, where we also obtain a particular solution describing 
the evolution a bremsstrahlung initial spectrum including a low-energy cutoff. Specific examples of the time-dependent energy 
spectra emitted by Comptonizing clouds are computed in § 5, and the implications of our results for the interpretation of 
X-ray data are discussed in § 6. Supplemental technical details of the mathematical approach are provided in Appendix A. 



2 TIME-DEPENDENT COMPTONIZATION 

In radiation-dominated, fully-ionized plasmas, photon creation and destruction are unable to establish local thermodynamic 
equilibrium, and the photons and electrons interact primarily via Compton scattering. In this section, we discuss the general 
equations governing the photon distribution as a function of time, space, and energy. The evolution equation for the photon 
energy distribution is solved in § 3. 



2.1 Transport equation 



Neglecting absorption, the temporal evolution of the distribution of photons in a nonrelativistic plasma cloud with electron 
number density n e and electron temperature T e is governed by the transport equation (Katz 1976; Payne 1980) 



dn 
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where n(e,f,t) is the photon occupation number (a dimensionless quantity), e is the photon energy, t is time, r is the 
spatial location within the plasma, a T is the Thomson cross section, c is the speed of light, m e is the electron mass, and 
k is Boltzmann's constant. The specific intensity, I v tx ergs cm -2 s _1 ster -1 Hz -1 , is related to the occupation number by 
I v = {2hv 3 1 'c 2 ) n, where h is Planck's constant and v — e/h is the photon frequency. The terms on the right-hand side of 
equation Q represent thermal Comptonization, spatial diffusion, and photon sources, respectively. We shall assume that /„, 
n, and the source term j are all isotropic. In this paper we will focus primarily on the process of thermal Comptonization 
occurring in homogeneous plasmas with steady properties. Effects related to bulk motion and temperature variations will be 
discussed in § 6. 

The occupation number n can be integrated with respect to e to obtain the total radiation number density, given by 



n r {r,i) 



^3 e n{e,r,t)de 



(2) 



The n 2 term in equation (0 describes stimulated scattering and is rarely important in astrophysical applications. When this 
term is neglected, the transport equation is rendered linear, and in the case of an isothermal plasma, we can obtain the 
equivalent form 
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where we have introduced the dimensionless photon energy 
e 



'«)= — • (4) 
By transforming the variable of integration in equation J2J from e to x, we can reexpress the photon number density as 



n(x, f, t) dx 



(5) 
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Since equation © is linear in n, it is sufficient to consider a source term j — j, that is localized in time, space, and energy, 
given by 



^8^1 (iS) S(t-t )5(r-r )5( X -xo) 



(6) 



which represents the injection of a single photon with dimensionless energy xo = eo/(kT e ) at time to and location fo. The 
solution to equation corresponding to the soure term j* is the Green's function, n G (x, xo, f, fo,t, to), describing the process 
of time-dependent thermal Comptonization in an isothermal plasma. 

Once the Green's function has been determined by solving equations © and JSJ in combination with a suitable boundary 
condition imposed at the surface of the cloud, the response to a general source term j can be obtained by performing the 
convolution 



n G (x, x , f, f ,t, t ) j(x ,t ,fo) x% dx dt d 3 f , 



(7) 



where the fo integration proceeds over the entire volume V of the cloud, regardless of its shape. Equation (|7J implicitly 
assumes that no radiation is present in the cloud at time t — 0, and that no radiation is incident on the cloud from the 
outside. As an alternative to solving for n G using equations © and @, we point out that n G is also the solution to the 
homogeneous equation 
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The homogeneity of equation ® makes it a more convenient starting point than equation © for the determination of the 
Green's function n G . In a homogeneous plasma with steady properties, both T e and n e are constants, and it is convenient to 
rewrite equation JSJ as 



dn G 
dy 



J_d_ 

x 2 dx 



i / dn„ 



+ V- 
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where we have introduced the dimensionless time 

kT e 



y(t) = a(t- to) 



(10) 



(11) 



This is the familiar Compton y-parameter, which generally must exceed unity in order for significant modification of the 
spectrum to occur (Rybicki & Lightman 1979). The constant a is the "Comptonization rate," and y — at the initial 
time t = to. In the case of a homogeneous plasma with constant density and temperature, the temporal evolution of the 
photon distribution depends only on the elapsed time since injection, t — to, and therefore we can write n G (x, xo,f, fo, t, to) = 
n G (x,xo,f,fo,y) without loss of generality. 

By operating on equation 1101 with 8Tv(kT e /hc) 3 x 2 dx, we can show that the radiation number density associated 
with the Green's function, 



v(r,fo,y) = 8ty (^—^j j x 2 n G (x,x ,f,fo,y)dx 
satisfies the homogeneous spatial diffusion equation 



= V- 



dy y3an e a T 



(12) 



(13) 



Note that the Comptonization term in equation 1101 vanishes upon integration, which is a manifestation of the fact that 
Compton scattering does not create or destroy photons. The initial condition for r\ can be obtained by combining equations © 
and 1121 1. which yields 



v(r,fo,y) 



y=0 



= 6(f- f ) ■ 



(14) 



Note that equations 113H and 114H must be supplemented by a suitable flux boundary condition imposed at the surface of the 
cloud in order to account for the escape of radiation from the plasma. 
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2.2 Separability 

We shall focus on the evolution of the occupation number Green's function n G subject to the initial condition given by 
equation The transport equation 1101 is separable in time, space, and energy when n e and T e are constants, as assumed 
here (Payne 1980). We employ separability by writing the occupation number n as the product 

n G (x,x ,r,r ,y) = ^- (x£r) nif,r ,y) f G {x,x ,y) , (15) 



•OV' -»"'>'<""- 87r \ kTe 

where f G (x,xo,y) is the Green's function for the photon energy distribution. Note that the energy and spatial dependences 
have been separated since the former is contained in f G and the latter is contained in 77. Substituting for n G in the transport 
equation 1101 using equation (1151 , and utilizing the fact that 77 satisfies equation 1131 , we can show that the energy distribution 
f G satisfies the Kompaneets (1957) Fokker-Planck equation, 



df G = j_d_ 

dy x 2 dx 
along with the initial condition 



(16) 



f G {x,x ,y) 



= x 8{x - xo) . (17) 

v=o 



The terms proportional to f G and df G /dx inside the parentheses on the right-hand side of equation 1161 express the effects 
of electron recoil and stochastic (second-order Fermi) photon energization, respectively. Due to our neglect of stimulated 
scattering, the asymptotic solution to equation 1161 obtained as y — > 00 is the Wien spectrum f G oc e~ x , rather than a Bose- 
Einstein distribution. The function f G (x,xo, y) represents the solution to the Kompaneets equation in an infinite, homogeneous 
medium. Spatial effects associated with the geometry and the finite size of the scattering cloud are introduced via the density 
distribution rj(f,fo,y). Note that n G satifies the initial condition given by equation @ as required. 

Once the spatial distribution has been obtained by solving equation 1131 for ri(f, fo, y), we must next determine the energy 
distribution f G (x,xo,y) in order to construct the complete solution for n G using equation 1151 . Although Payne (1980) and 
Rybicki & Lightman (1979) state that equation 1161 must be solved numerically in general, we shall demonstrate below 
that an analytical solution for the Green's function can in fact be obtained in the form of a real integral. By operating on 
equation 116^ with J °° x 2 dx, we can establish that f G has the convenient normalization 

x 2 f G (x, xq, y) dx = constant = 1 , (18) 



i) 



where the final result follows from the initial condition (eq. 1171 '). Note that this normalization is maintained for all values of 
y, which reflects the fact that Compton scattering conserves photons. We shall seek to solve equation 1161 using Laplace 
transformation in § 3. Once the solution for the Green's function f G (x,xo,y) is known, the particular solution for the 
distribution function f(x,y) corresponding to an arbitrary initial spectrum fo(x) can be found using the integral formula 



f(x,y) = x fo{x )f G (x,x ,y)dxo. (19) 
Jo 

By operating on equation (11911 with J"°° x 2 dx, we can establish that 

poo poo poo 

h(y)= / x 2 x 2 l f (x )f G (x,x ,y)dx dx = / x 2 f {x ) dx = constant , (20) 



where the final result is obtained by reversing the order of integration and applying equation 1181 . The constancy of I2 follows 
from the fact that Comptonization does not create or destroy photons, and we therefore refer to I2 as the "number moment" 
of the particular solution f(x, y). Equation 1201 will provide a useful check when we consider the properties of the particular 
solutions obtained in 8 4. 



3 SOLUTION FOR THE GREEN'S FUNCTION 

Our analytical approach to the determination of the Green's function energy distribution f G (x, xo,y) will be based on Laplace 
transformation. For the sake of brevity, only the primary steps are discussed here. Readers interested in the technical details 
are referred to Appendix A. 

3.1 Laplace transformation 

Laplace transformation of the Kompaneets partial differential equation 1161 with respect to y yields the ordinary differential 
equation 
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S L - X 2 5(x - X ) = \ 4~ 

ar ax 



(21) 

Jb~ LLJb L \ UX / J 

where 

L(x,x ,s)= e~ sy f G (x,x ,y)dy (22) 



denotes the Laplace transform of f G (x, xo, y), and we have also used equation l|17fl . We show in Appendix A.l that the solution 
for the transform L(x,xq,s) is given by 

L(x, x ,s) = r r (M ~ 3 2 /2) xo 2 x~ 2 e^' 2 (23) 
1 W {M 2 , fl (x )W 2 , fi {x) , x>x , 

where the quantity fi is a function of the transform variable s, defined by 

/ 9\ 1/2 

Ks) - (24) 

and M2, fi(x) and W2, ^(a;) denote the Whittaker functions (Abramowitz & Stegun 1970). Equation l!2.St can be rewritten in 
the equivalent form 



L(x,x ,s) = ^ + ^ 2) xo 2 x~ 2 e ^°- x)/2 M 2 , M (x min ) W 2 , M (x max ) 



(25) 



where 

imin = min(x,a;o) , x max = max(i, x ) ■ (26) 
The function W2, n(x) is defined in terms of M2,p,(x) by equation (13.1.34) of Abramowitz & Stegun (1970), which gives 

W^(x) . r( ^ /2) M^(x) + f ^- M ^) , (27) 

and the function M2, i_i(x) can be evaluated using the series expansion 



M 2 ,^x) = e- x/2 x» +1/2 



M-3/2 (a. -3/2)0* -1/2) 



(28) 



l + 2^i (l + 2/Lt)(2 + 2/x) 2 

The solution for Z/(x,o;o,s) in equation 1251 has been obtained by utilizing continuity and derivative jump conditions based 
on the differential equation 121H . along with various identities satisfied by the Whittaker functions. Further details of the 
derivation are provided in Appendix A.l. 

3.2 Integral expression for the Green's function 

To obtain the solution for the Green's function / G (x, xo, y), we must perform the inverse Laplace transformation of L(x, xo, s) 
using the complex Mellin inversion integral (Butkov 1968), 



1 



f G (x,x ,y) = I e sy L(x,x ,s)ds , (29) 

J ~y — ioo 

where the constant 7 is selected so that the line Re s = 7 lies to the right of any singularities in the integrand. In this problem, 
singularities occur where the quantity fi — 3/2 is zero or a negative integer, resulting in the divergence of T(n — 3/2). The 
definition of /it (eq. [24] ) therefore implies that simple poles are located at s = and s = —2. The corresponding values for 
H are fx = 3/2 and /i = 1/2, respectively. From the locations of the singularities, it follows that we must require 7 > for 
convergence of the integral in equation l|29^ . 

The simple poles at s = and s = —2 both lie to the right of the branch point for the square root function, which is 
located at s — —9/4. Hence they are contained within the closed integration contour C indicated in Figure 1. We can therefore 
use the residue theorem to write 

2 

e sy L(x, xq, s) ds — 2ni Res(s„) , (30) 



c 



n = l 



where the left-hand side denotes the integral around the contour C and Res(s,j,) is the residue associated with the simple pole 
located at s = s n , with si = and s 2 = —2. Note that the integration contour must avoid the branch cut of the square root 
function, extending from s — —9/4 to s = —00. Asymptotic analysis indicates that the contributions to the integral along the 
large arcs MN and QR vanish in the limit n — > 00 (see Fig. 1). Likewise, the integration along the small arc OP vanishes in 
the limit r 2 — * 0. Our solution for the energy distribution / G therefore reduces to 
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Figure 1. Closed integration contour C for the Laplace inversion integral on the left-hand side of equation 1301 . The branch point of 
the complex square root function is located at s = —9/4, in the center of the small arc OP, which has radius T2- Simple poles are located 
at s = and s = —2; these singularities are contained within the contour C provided 7 > 0. The branch cut (dot-dashed line) extends 
from the branch point to s = —00. See the discussion in the text. 



1 f° 1 f Q 

f G (x, x , y) — —■— / e sy L(x, x , s) ds - ■— / e 3y L(x, x , s) ds + > Res(s n ) , (31) 

2m J N 2m Jp fri 

in the limit fj — > 00, 7-2 — > 0. Equation 131^ expresses the Laplace inversion in terms of two residues and two integrals, one 
above the branch cut and one below it. 

The residues appearing in equation 1311 are evaluated in Appendix A. 2, where we find that (see eqs. IA18I ) 

n=l 

The procedure for treating the integrations above and below the branch cut is discussed in Appendix A. 3. This involves 
the utilization of symmetry relations and other identities satisfied by the Whittaker functions. The final result given by 
equation 1A29|I is 

32 -9y/4 -2 -2 (x„-x)/2 f°° -u 2 y U Smh(7Tu) 



f G (x,x ,y) = - e-Wxo'x-'e^-*" Jo , (1 + 4?t 2 )(9 + 4m 2) 

X W 2 , l u{x )W2, l u{x)du + + — ~r~~~ ^ - ^ ^ " ^ . (33) 

2 2 Xo X 

This solution for the Green's function is one of the main results of the paper. It represents the fundamental "kernel" for the 
process of thermal Comptonization in an infinite homogeneous medium, from which particular solutions for any distribution 
of photon sources can be developed by quadrature. Spatial effects in a finite medium can be incorprated using the separation 
solution given by equation I15H . In Appendix A. 4, we provide the series expansions needed to evaluate the Whittaker functions 
in equation l|33^ . Furthermore, a self-contained FORTRAN code that evaluates the Green's function by performing the 
integration in equation 133^ is available from the author upon request. In the remainder of the paper, we will analyze the 
properties of the Green's function, and demonstrate that all of the previously known steady-state and time-dependent solutions 
can be reproduced using it. Hence equation 133H represents a significant generalization of the earlier analytical results in the 
theory of time-dependent thermal Comptonization. 
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x 

Figure 2. Analytical solution for the Green's function x 2 f G (x, xo,y) (eq. 1331 : solid lines) plotted as a function of the dimcnsionlcss 
photon energy x for the indicated values of the dimensionless time y. The area under each curve is equal to unity, and the initial photon 
energy is given by (a) xq = 0.1; (b) xq = 1; (c) xq = 5; and (d) xo = 10. Included for comparison are the numerical solutions (dashed 
lines) obtained by integrating numerically the Kompaneets equation using a Gaussian initial condition, as explained in the text. The two 
sets of curves are essentially identical. As y increases, x 2 f G (x, xo,y) approaches the Wien spectrum (1/2) x 2 e~ x (filled circles) in each 
case as expected. 

3.3 Approach to Wien equilibrium 

In the time-dependent, isothermal Comptonization problem treated here, the Green's function must approach the Wien 
equilibrium spectrum, i.e., 



lim f G (x,xo,y) = f w (x) =7]^ 

y — >oo Z 



(34) 



where the factor of 1/2 on the right-hand side is required in order to ensure that f w satisfies the normalization condition 



x /w ( x ) dx ■■ 



(35) 



in compliance with equation 1181 . Our solution for f G (x,Xo,y) expressed by equation 133H clearly satisfies equation 134H . 
and therefore the Green's function exhibits the correct asymptotic behavior for large values of y. In Figure 2 we plot f G as 
a function of the dimensionless photon energy x and the dimensionless time y for four values of the initial photon energy, 
xo — 0.1, xo = 1, xo — 5, and xo = 10. Also included for comparison is the Wien equilibrium spectrum, f w (x) = (1/2) e~ x . 
Note that in each case, the Green's function approaches the Wien shape for large values of y as predicted. The spectrum 
depicted in Figure 2 describes the Comptonization of an intrinsically narrow line feature. Applications to the broad Fe Ka 
lines observed in AGNs are discussed in S 6. 



3.4 Comparison with numerical simulations 

The acid test for our Green's function solution is provided by comparing it with the photon energy distribution obtained by 
solving numerically the Kompaneets partial differential equation 11611 . In addition to specifying the initial condition for the 
photon distribution, the numerical solution procedure also requires the imposition of boundary conditions for f G at large and 
small values of x, for all values of y. Obviously, this information is not available a priori, and therefore we are forced to push 
the boundaries of the computational domain as far as possible away from the region of interest, so as to avoid contaminating 
the solution with imprecise boundary data. This requirement markedly reduces the efficiency of the numerical approach, since 
much of the computer time is spent on calculations outside the region of interest. We have performed numerical simulations 
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using the IMSL partial differential equation solving subroutine DMOLCH. The initial condition is specified as a narrow 
Gaussian feature, with selected values for the mean energy (a;) and standard deviation a. This initial condition approximates 
the 5-function distribution given by equation 1171 when xq = (x) and <j/{x) <C 1. The analytical and numerical solutions for 
f G (x,xo,y) are compared in Figure 2 for four different scenarios. In Figure 2a, we set (x) = 0.1 and a = 0.01; in Figure 2b 
we set (x) — 1.0 and a — 0.01; in Figure 2c we set (x) = 5.0 and a = 0.05; and in Figure 2d we set (a;) = 10 and a = 0.1. The 
results are plotted for several values of y. Note that the numerical and analytical solutions are virtually indistinguishable, 
confirming that our closed-form expression (eq. 1331 ) is in fact an exact representation of the Green's function. 



4 PROPERTIES OF THE GREEN'S FUNCTION 

The solution we have obtained in § 3.2 for the Green's function / G (a;, xo,y) given by equation l|33|l has a number of interesting 
properties that are further explored in this section. 



4.1 Moments of the Green's function 

We can obtain additional insight into the behavior of the Green's function by focusing on the variation of the "power moments," 
In (y), defined by 

poo 

In(y)= / x n f G (x,x ,y)dx . (36) 



Although the results obtained below for In(y) are valid for general real values of n, in our specific applications we shall focus 
on integral values, since these correspond to cases of special physical interest such as the "number moment" I§ introduced 
in equation l|20|l . and the "energy moment" I§ discussed below. In the case of the Green's function, the initial values of the 
power moments at y — are obtained by substituting equation 1171 into equation 136H . which yields 

7°(0) = xr 2 ■ (37) 
We can determine the subsequent evolution of I„{y) by operating on equation l|33[l with x n dx. The result obtained is 

T G M _ 32 9 „/4 -2 Xo /2 ^ -u'y " sinh(^) 

l n {y)-_e x e / e (i + 4w 2 )(9 + 4w2) ^,»(x j 



;T JO 

-x/2n-2 u/ ,„ w _ j„. , r(n+l 



' W 2 , iu{x) dx du + r(n n + 1} + (2 - n) T(n) e~ 2y (— - I) , (38) 

where we have interchanged the order of integration over u and x. The integration with respect to x converges provided 
n > 1/2, in which case we can use equation (7.621.11) from Gradshteyn & Ryzhik (1980) to find that 

rG^A _ 32e- 9y/4 e xo/2 f°° - u * y u sinh(Tru) 



™= ;,gr ( n- 2) I e "" a (TT^#Tl^^^ )r(n - 1/2 + ^ 



X T(n - 1/2 - iu) du + F(n+1) + (2 - n) T(n) e~ 2y f— - -) . (39) 

2 V Xo 2 / 

We shall examine the behavior of this general expression for several cases of special interest below. 



4.2 Variation of the number and energy moments 

We have established via integration of the Kompaneets equation I16t that / G (a;, xo, y) must satisfy the normalization condition 
(cf. eq. ED) 

l2{y)= / x 2 f G (x,x ,y)dx = l (40) 



for all values of y, in agreement with the initial condition for I§{y) expressed by equation 1371 . As pointed out in § 2.2, this 
behavior is a consequence of the fact that Compton scattering conserves the photon number density. We can use equation I39t 
to confirm explicitly that this condition is actually satisfied by our solution for the Green's function. When n = 2, the integral 
in equation l|39^ vanishes due to the factor T(n — 2) in the denominator, and we obtain 

i2(y) = i , (41) 

in satisfaction of equation 1401 . We have therefore established that our Green's function solution conserves photons and is 
properly normalized. Next we focus on the variation of the "energy moment," 
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Figure 3. Mean photon energy x for the case of a monoenergetic initial spectrum (eq. 1461 ). plotted as a function of the dimensionless 
time y for the indicated values of the initial photon energy Xq. Note that in each case, x — > 3 as y — > oo because the spectrum always 
equilibrates to the Wien form. The photons experience continual heating, on average, when xq < 2, and continual cooling when xq > 4. 
In the intermediate case with 2 < xq < 4, a period of initial heating is followed by cooling. This behavior is a manifestation of the 
underlying variation of the spectral shape, as explained in the text. 



i^(y) 



x f G (x,x ,y)dx , 



which is related to the radiation energy density, U r , via 

(fcT e ) 4 



U r (r,t) = J — e 3 n{e,f,t)de 



8vr- 



(ch) s 



if GO 



Note that 1$ (y) is also equal to the mean photon energy, x(y), since 
-t \ $™ ^ j G {x,x ,y)dx 

x (y) = r oo : . — - = h (y) , 



(42) 



(43) 



(44) 



J™ x 2 f G ( x > x o,y)dx 

where the final result follows from equation I4H . The variation of the mean photon energy was also considered by Kompaneets 
(1957), and we shall compare our result for x(y) with his. Setting n = 3 in equation 1391 and using the identity 

tt(1 + 4u 2 )(9 + 4« 2 ) 



r(5/2 + iu) T(5/2 - iu) 

16 cosh(7rtt) 

we obtain for the variation of the mean photon energy 



Ky) = ^3 G (2/) = 2: 



-9y/4 e x /2 



W 2 , 



Axo) u tanh(7rit) du + 3 — 2e 2y ( — —~\ 

\xo 2 / 



(45) 



(46) 



This result agrees with equation (48) from Kompaneets (1957), aside from the fact that his expression for x(y) contains an 
erroneous factor of e~ x °^ 2 in place of the correct factor e x °^ 2 appearing in equation 146|l . In Figure 3 we plot the variation 
of x(y) as a function of y for several values of the initial photon energy xq. Note that in all cases, we find that x(0) = xq 
initially, as required. For large values of y, we find that x — > 3, as expected based on the fact that the photon distribution 
should approach the Wien spectrum f w (x) oc e~ x . However, the evolution of x is not necessarily monotonic. For example, 
when xo = 3, the mean energy x increases until y ~ 0.25, and then it decreases, asymptotically approaching the Wien value 
of 3 as y — * oo. This behavior can be understood by examining the associated evolution of the inverse-Compton temperature, 
which we investigate in § 4.3. 



4.3 Variation of the inverse-Compton temperature 

By operating on the Kompaneets equation 116H with x n dx and integrating by parts twice, we can show that the power 
moments In(y) satisfy the "differential recurrence relation" 
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^- = (n-2)[(n+l)/°(j/)-/° +1 (j,)] . (47) 
In particular, setting n = 3 yields the conservation equation for the energy moment lf(y) = x{y), 



, T 10 {y) 



(48) 



(49) 



where T IC (y) denotes the inverse-Compton temperature of the radiation field, defined by 

T IC (y) _ l l?(y) 

T. ~ 4 If(y) ' 
and 

p oo 

I?(y)= / x 4 f G (x,x ,y)dx . (50) 



Equation Ij48^ describes the conservation of total energy for a group of photons experiencing thermal Comptonization in a 
plasma with a Maxwellian electron distribution. When T IC < T e , the photons gain energy (x increases), and x decreases 
when T IC > T e . The inverse-Compton temperature varies in response to changes in the shape of the spectrum, and it does 
not necessarily track the behavior of x. No net energy is exchanged between the radiation and the matter when T IC — T e , 
although subsequent evolution of the spectral shape can cause T IC to move away from T e . In the limit y — > oo, we find that 
T IC — > T e regardless of the initial energy of the photons because the spectrum asymptotically approaches the Wien form 
/ w oc e~ x . 

In the case of a monoenergetic initial spectrum, we can use equations 1371 1 and 1491 to conclude that the initial temperature 
ratio is given by 



T IC (0) _ x 



(51) 



T e 4 

Hence the photons will initially gain energy if xo < 4, and otherwise they will lose energy. The subsequent time evolution of 
T ic (y) can be determined by computing the variation of lf(y), which is accomplished by setting n = 3 in equation 14711 and 
substituting for lf(y) using equation 1461 . The result obtained is 

p-Sy/4 p xo/2 r°° 2 / 1 1 \ 

!?(y) = 7T-2 / e ~" " W», iu{x ) u tanh(Tm) (25 + 4 M 2 ) du + 12 - 12 e - 2y ( . (52) 

2x o Jo Vx ° 21 

It is interesting to note that the same result can also be obtained by setting n — 4 in equation 1391 . which provides a useful 

check on the self-consistency of our formalism. Taken together, equations l|46|l . I|49|l . and ilo2li provide a new, closed-form 

solution describing the variation of the inverse-Compton temperature for an initially monoenergetic photon distribution. This 

result has not appeared previously in the literature, and it can be used to obtain additional physical insight into the energetics 

of thermal Comptonization, as discussed below. Our results for 1^ and 7^ confirm that T ic /T e — * 1 as y — > oo, since /.p — + 12 

and ^3* — > 3. This behavior is consistent with our expectation that the spectrum should equilibrate to the Wien form. 

In Figure 4 we plot the results obtained for T IC (y) by combining equations 1461 , 1491 , and 1521 for several different values 

of the initial photon energy xo- There are a number of interesting features in the plots. First, note that for the cases with 

xq < 4, the photons are initially cooler than the electrons (see eq. [51)1 . and this explains the initial increase in the mean 

energy x displayed in Figure 3. Conversely, when xo > 4, energy is initially transferred from the photons to the electrons, 

and therefore x decreases. Note, however, that the variations of x(y) and T IC (y) are not necessarily monotonic. For example, 

the case with xo — 3 displays initial heating, and x reaches a peak at y ~ 0.25. Beyond this point, the inverse-Compton 

temperature exceeds T e , and the photons begin to lose energy. This is a consequence of the fact that x must eventually 

decrease to the Wien value of 3 as y — > oo, and this decrease can only happen if T 1C > T e . Ultimately, T 1C approaches T e 

and equilibrium is achieved. Hence the behaviors of T 1C and x are fundamentally driven by the underlying evolution of the 

radiation spectrum. 

For initial photon energies in the range 4 < xq < 6, equation 1511 indicates that T IC exceeds T e initially, and therefore 
the photons must lose energy to the electrons. However, despite the fact that x is initially decreasing, we observe that T IC 
actually proceeds to increase. This apparently paradoxical behavior stems from the variation of the spectral shape. We can 
obtain some useful insight into this phenomenon by computing the initial value of the temperature derivative, dT lc /dy, in 
the case of a monoenergetic initial spectrum. By differentiating equation 1491 with respect to y and using equation 1471 to 
evaluate the derivatives of the moments, we find that for general values of y, 

±^_*&_1# + 1(#\\ (53) 



dy T e 2 If 2 If 4 \lf 

The initial value of the temperature derivative at y = can be obtained by using equation 1371 to substitute for the moments 
which yields 
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Figure 4. Inverse-Compton temperature ratio T ic (y)/T e = (1/4) {y)/I§ {y) plotted as a function of y for the case of a monoenergetic 
initial spectrum with the indicated value of the initial photon energy xo . The functions 1^ (y) and 7jp ( y) are evaluated using equations 1461 
and 1521 . respectively. The inverse-Compton temperature increases initially if xq < 6. 



d T w (y) 



dy T e 



y = 



X ° Id \ 

— (6-a;o) 



(54) 



This result clearly demonstrates that the inverse-Compton temperature T IC initially increases if and only if xq < 6. On the 
other hand, according to equations l|48|l and l|51|l . the mean photon energy x initially increases if and only if xq < 4. Hence 
there is an intermediate regime, 4 < xo < 6, within which x initially decreases, while T 1C simultaneously increases. This 
illustrates the fact that x and T 1C are not directly connected, but each is determined by the underlying evolution of the 
spectral shape, which is governed by the Kompaneets equation. 



4.4 Particular solution for Wien initial spectrum 

The Green's function solution given by equation can be used to obtain the particular solution f{x, y) corresponding to 
any desired initial spectrum fo(x). As discussed in § 2.2, the particular solution is given by the convolution 

/•oc 

f(x,y)= xlfo(x )f G (x,x ,y)dx . (55) 
Jo 

The Wien initial spectrum, 

fo(x) = e~* , (56) 

is a case of special interest, since it represents the equilibrium solution to the Kompaneets equation 1 Kit . Substituting 
equation l|56|l into equation l)55|l and evaluating f G (x,xo,y) using equation we obtain after interchanging the order of 

integration over xo and u 



}{x,y) 



+ ^ e - 9y/4 x- 2 e- x/2 



u sinh(7ra) TTr . . f°° _,„«,„ , . , 



(57) 



The integration over xo in equation I57H can be performed using equation (7.621.11) from Gradshteyn & Ryzhik (1980), which 
yields zero. The particular solution in the case of the Wien initial spectrum therefore reduces to 



f(x,V) = e x , 

which confirms that the Wien spectrum remains unaffected by isothermal Comptonization, as expected. 



(58) 



4.5 Particular solution for bremsstrahlung initial spectrum 

The reprocessing of optically thin bremsstrahlung emission in clouds of hot electrons is thought to play a major role in the 
formation of the X-ray spectra observed during intense AGN flares (Lightman, Giacconi, & Tananbaum 1978). Becker & 
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Begelman (1986) and Chapline & Stevens (1973) studied this process by analyzing the evolution ol an initial spectrum of the 
form 

/<,(*) =aT 8 e— . (59) 
The exact particular solution to this problem is given by Becker & Begelman as 

f BS (x,y) ee x- 3 e-{l + f(l- e - 2H ) x + | [1 - e-^ (2y + 1)] x 2 + | [y - 1 + e" 2 ^ + 1)] x 3 | . (60) 

This solution demonstrates the distortion produced by electron scattering, and the development of a thermal peak in the 
spectrum. However, due to the fact that the initial spectrum contains an infinite number of low-energy photons in this case, 
it is not possible to achieve full equilibration to a Wien distribution even in the limit y — -> oo. This unphysical behavior is 
an artifact of the utilization of equation l|59[l to describe the initial spectrum all the way down to zero energy. In an actual 
astrophysical situation, self-absorption effectively introduces a low-energy cutoff in the initial spectrum, as discussed below. 

During intense AGN flares, soft seed photons produced in a relatively dense source region via bremsstrahlung are thought 
to be reprocessed in a corona containing hot electrons, in which scattering dominates over photon creation and destruction 
(Sunyaev & Titarchuk 1980; Nandra 2001). However, for photon energies below some critical energy x*, the mean free path 
for free-free absorption in the source region becomes smaller than the size of the source. If the source region is homogeneous, 
then we can use the usual formulas for thermal bremsstrahlung to show that (Rybicki & Lightman 1979) 

r 3 U 

6-6^7, (61) 



1 - e~ x * ' aT 4 

where T and U r denote the gas temperature and radiation energy density in the source region, respectively, and aT 4 is the 
blackbody energy density. In the physical applications of interest here, U r <C aT 4 , and therefore we can rewrite equation Hill 

as 

'T er x ' 



2.6 



(f ) « 1 , (62) 

where T e g = (Ur/a) 1 ^ 4 is the effective temperature of the radiation in the source region. Below the critical energy x*, the 
initial spectrum given by equation I60H transitions into a Planck distribution at temperature T due to self-absorption. Since 
the number of photons per unit frequency range contained in the Planck distribution is proportional to the frequency v as 
v — > 0, this effectively introduces a cutoff in the initial spectrum at energy x — x*. The existence of this low-energy cutoff 
has a profound effect on the spectral evolution in the scattering region because the number of photons contained in the initial 
spectrum is now finite. 

With the availability of our analytical expression for the Green's function given by equation l|33[l . we are in a position 
to improve the situation by treating an optically thin bremsstrahlung initial spectrum that includes a low-energy cutoff. 
Specifically, we shall analyze the evolution of the "modified" bremsstrahlung initial spectrum given by 

ill: ; (63) 

where x* < 1 denotes the low-energy cutoff, which approximates the effect of self-absorption. The number moment associated 
with this initial spectrum is (see eq. 1201 ') 

/>oo 

h= / x 2 /»(# = r(0,i,) , (64) 
Jo 

where T(a,z) denotes the incomplete gamma function (Abramowitz & Stegun 1970). Note that the modified bremsstrahlung 
initial spectrum contains a finite number of photons. Combining equations 1331 . I55II . and 1631 1. and interchanging the order 
of integration over xo and u, we obtain 

f(x,y)=™e- 9 ^x- 2 e-^ H e -»'» WiMx) [°° x ^ e^W^o) dx du + S, (65) 

where 

S=\e-'T(0,x.)+e-"-*' Q - |) [2T(-1, x.) - T(0, x.)] . (66) 

The integration over xo in equation 1651 can be worked out by employing equation (13.1.33) from Abramowitz & Stegun 
(1970) and equation (3.2.12) from Slater (1960) and integrating by parts three times. After some simplification, we obtain for 
the particular solution 

ti \ 32 -9y/i -2 -2 -ix+x,)/2 f 00 -u 2 y u sinh(nu) 



tt J (1 +4w 2 )(9 + 4m 2 ) 

x [W r i,,u(x.)-3Wo,iu(x.) + 6W_,, iu (x.)-6W-2,iu(x.)] du + S. (67) 
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Figure 5. Time-dependent photon energy distribution x 2 f(x, y) resulting from the Comptonization of the modified bremsstrahlung 
initial spectrum (eq. 1671 : solid lines) plotted as a function of x for the indicated values of the dimensionless time y. Also included for 
comparison are the corresponding initial spectra (eq. 1631 : dotted lines) and the equilibrium Wien spectra (filled circles). The values for 



the low-energy cutoff in the initial spectrum are (a) x* 



10" 



(b) x, 



10" 



(c) X* 



10" 



(d) x* = 0. In panel (d), we have used 



equation I6UI to evaluate as /, and the Wien spectrum has been omitted since it does not exist in this case. The solutions in panels (a) 
- (c) show proper equilibration to the Wien spectrum as y increases, whereas the distribution in panel (d) does not due to the infinite 
number of photons in the initial spectrum. 



This result describes the time-dependent Comptonization of an optically thin bremsstrahung initial spectrum with a low- 
energy cutoff at x = x* and a finite number of photons. Note that as y — > oo, the integral vanishes and S — > r(0, x*) e~ x /2, 
and therefore we find that 



lim f(x,y) 



T(0,x«) 



(68) 



where I2 = T(0, a;*) according to equation 1641 1. Hence our solution clearly equilibrates asymptotically to the Wien spectrum 
/w ( x ) = (^2/2) e ~ x 1 which satisfies the normalization requirement f 00 x 2 f w (x) dx = J2. Series expansions useful for evaluating 
the Whittaker functions appearing in equation 16711 are provided in Appendix A. 4, and a self-contained FORTRAN code that 
performs the integration to determine f(x,y) is available from the author upon request. 

The particular solution given by equation 1671 represents a generalization of the x* = result f BB (x,y) (eq. I6UI ). and 
we find that 



lim f(x, y) = / BB (x, y) 



(69) 



as required. In Figure 5 we use equation 1671 to plot f(x,y) for several values of the low-energy cutoff x*. We also depict the 
equilibrium Wien spectrum f w (x) = (I2/2) e~ x , with I2 ~ F(0,a;,). As expected, the solutions with nonzero values of a;* are 
able to fully equilibrate to the Wien spectrum because the number of low-energy photons is finite in these cases. Equilibration 
occurs more rapidly as x, is increased. Also included for comparison is the solution / BB (x, y) corresponding to x* = 0, with 
an infinite number of photons (eq. 1601 ). In this case, the spectrum is unable to equilibrate. Our new solution for f(x,y), 
given by equation 1671 . therefore displays a more reasonable physical behavior, and consequently it provides a rigorous basis 
for models of the X-ray spectral evolution during AGN flares. In § 5 we discuss the results obtained when equation 1671 is 
coupled with a spatial transport model. 



5 MODELS INCLUDING PHOTON ESCAPE 

The Green's function solution obtained in § 3 describes the evolution of the photon energy distribution f G (x,xo,y), but it 
does not address spatial issues such as the distribution of photons within the scattering cloud. In order to obtain the Green's 
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function for the occupation number, n G (x, xo, r, fo, y), we must solve the spatial diffusion equation 113H for the radiation 
number density r/(f, fo, y), and combine this information with / G (x, xq, y) using equation 115H . However, in many applications 
the detailed spatial distribution is not needed, and in such cases one may simply integrate over the volume of the cloud 
to describe the radiation spectrum as a function of energy and time only. In this section we take the latter approach and 
incorporate spatial effects by using a simple escape-probability formalism that serves to illustrate the general concept. More 
complex models with realistic geometries can easily be accommodated within our framework, provided the electron distribution 
in the cloud is homogeneous, as we have assumed. 



5.1 Escape time distribution 

Let us suppose that iVo photons are injected into the cloud with an arbitrary spatial distribution at time t = to- Then 
the probability that a randomly-selected photon still remains inside the cloud at time t > to is given by the escape time 
distribution 

"<«>-¥■ <™> 

where 

N r (t) = j> n r (f,t)d 3 f (71) 

denotes the number of photons contained inside the cloud at time t, which is computed by integrating the radiation number 
density n r over the volume V of the cloud. It follows that P(t) decreases monotonically as a function of time from unity at 
t — to to zero as t — > oo due to the escape of photons from the cloud. 

In a homogeneous plasma, the energy and spatial distributions evolve in a decoupled manner as demonstrated in § 2.2. 
The probability that a photon randomly selected from the iVo photons injected at time t = to will subsequently escape in the 
time interval between t and t + dt can be expressed as 

„, . dt 1 dN r , . . 

where the escape probability distribution P(t) is given by equation (ZD) ; and 
1 dPN 



fo 



*»--(?¥) < 73 > 

represents the mean escape time for photons remaining in the cloud at time t. Note that P(t) dt/t csc is also equal to the 
fraction of the initial photons escaping between times t and t + dt. It is important to emphasize that, in general, the mean 
escape time f eS c varies as a function of time t due to the evolving nature of the photon distribution inside the cloud, which 
is a result of spatial diffusion. This variation can occur even when the electron distribution in the cloud is homogeneous, as 
assumed here. By virtue of equation l|72^ . P(t) clearly satisfies the normalization requirement 

P(t) -* = 1 , (74) 

£esc 
*0 

since N r (to) = No. Once the radiation number density distribution n r (f,t) has been determined by solving the appropriate 
diffusion equation, the escape probability P(t) can be evaluated using equations I7U1 and 17111 . The result obtained for n r (f,t) 
will naturally depend on the shape of the cloud and also on the initial spatial distribution of the photons at time t = to. 
Consequently, these factors will also influence the escape time distribution P(t). Interested readers can find further details in 
Sunyaev & Titarchuk (1980). 

Let us suppose that P(t) has been computed in a given situation, and that the particular solution for the photon energy 
distribution f(x,y) has also been determined by utilizing equation H19I . In this case, the spectrum of the radiation remaining 
in the plasma at time t can be obtained by forming the product of / with the probability distribution P. Specifically, we find 
that the number of photons with dimensionless energy between x and x + dx still residing within the cloud at time t is given 
by 

N x (x,t)dx = P(t) x 1 f(x,y)dx , (75) 

where y(t) = (t — to)n e a T c(kT e /m e c 2 ) (see eq. Hil l. Integration of N x with respect to energy yields the total number of 
photons in the cloud, 

N r (t) = / N x {x,t)dx . (76) 







The number of photons with energy between x and x + dx escaping from the cloud in the time interval between t and t + dt 
is then given by 
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N x {x,t)dtdx = t cs l N x (x,t)dtdx = t~ 1 c P(t) dt x 2 f(x,y) dx . (77) 

In § 4 we obtained particular solutions f(x, y) corresponding to several specific initial spectra. The associated results for the 
escaping radiation spectrum N x are discussed below. 



5.2 Escaping spectrum for monoenergetic initial condition 

In order to demonstrate the effects of spatial transport without undue complexity, we shall assume here that t csc — constant. 
This implies that the spatial distribution of the photon sources is proportional to the first eigenfunction of the diffusion 
operator, as discussed by Sunyaev & Titarchuk (1980). The specific result obtained for the initial density distribution in the 
spherical, homogeneous case with t csc — constant is 



n r (r,t) 



No A 2 sm(\r/£) 



sm{\R/£) - (XR/£) cos(A7?/£) A-kP r 



(78) 



where No is the total number of photons in the initial distribution, R is the radius of the cloud, £ = (n e a T ) 1 is the (constant) 
mean free path for electron scattering, and the eigenvalue A is the smallest positive root of the equation 

tan (-) = (-) I 1 --) • (79) 

The quantity a introduced in equation l|790 is a constant of order unity that depends on the precise manner in which the 
flux boundary condition is specified at the surface of the cloud using the Eddington approximation. For example, Rybicki & 
Lightman (1979) set a = 3 1 / 2 in their equation (1.124), whereas Sunyaev & Titarchuk set a = 3/2 in their equation (A. 3). 

As t increases from the initial time to, photons escape from the cloud due to spatial diffusion. The resulting time-dependent 
solution for the radiation number density is given by 

n r (r,t) =n r {r,t )e- {t - to)/t ^ , (80) 
where n r (r,to) is the initial distribution [cq. I|78|l ] and 
31 

Use = T-5— = constant (81) 
A c 

denotes the mean escape time. The probability that a photon injected at time to still remains within the cloud at time t is 
therefore given by 

P(t) = e-C-'oJ/w _ ( g2 ) 

We can now use equation 175H to express the number of photons with dimensionless energy between x and x + dx residing 
within the cloud at time t as 

N x (x,t)dx = e -(*- 4 °)/w x 2 f( Xy y)dx , (83) 
where y(t) = y(t- t )/t csc , and 

kT e 2 kT e 

y = t esc n e a T c e — = Max(r, r ) e — (84) 



represents the mean value of the Compton {/-parameter experienced by photons before escaping from a cloud with scattering 
optical thickness r = Rjl (Rybicki & Lightman 1979). When y^l, significant Comptonization occurs before most of the 
radiation has escaped. In the limit of small y, little change in spectral shape occurs during the transient. The number of 
photons with energy between x and x + dx escaping from the cloud in the time interval between t and t + dt can be written 
as (see eq. 1771 ) 

N x (x,t)dtdx = t~l e - ( *- to)/w dt x 2 f(x,y)dx . (85) 

In Figure 6, we plot the escaping photon distribution N x (x,t) evaluated using equation I85|l . with f(x,y) set equal to the 
Green's function solution given by equation 133H . Results are presented for two values of the initial energy xo and two values 
of the mean Compton parameter y. The curves depicted in Figure 6 represent the spectra that would be observed outside 
the cloud as the initial burst of monochromatic (e.g., iron line) radiation is reprocessed by electron scattering, and gradually 
escapes from the cloud. For small values of t — to, the escaping spectrum is narrow and centered on the initial energy xo since 
little scattering has occurred. At intermediate times, the spectrum is broadened and the detailed shape depends on the value 
of the initial energy xo- The photon distribution is characterized by spectral hardening for xo<3 and by spectral softening 
for xo ^ 3. For large values of t — to, the spectrum has a Wien shape with a low normalization because most of the radiation 
has escaped. Equilibration to the Wien spectrum occurs more rapidly as xo and/or y are increased. 
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Figure 6. Time-dependent escaping photon distribution N x (x,t) (eq. 1851 : solid lines) plotted in units of ^ c for the case of a spectrum 
that is initially monoenergetic at time t = to- The elapsed time t — to in units of t CBC is indicated for each curve, and the values of the 
initial photon energy xo and the mean Compton parameter y (eq. 1841 ) are given by (a) xq = 0.1, y = 0.2; (b) xq = 0.1, y = 1; (c) 
xq = 10, y = 0.2; and (d) xq = 10, y = 1. The normalization of the curves decreases over time due to the loss of radiation from the 
scattering cloud, as described by equation 1821 . Included in each case is the Wien spectrum (filled circles) evaluated at time t max , where 
tmax is the maximum time for the sequence of curves in the panel. 

5.3 Escaping spectrum for bremsstrahlung initial condition 

It is also interesting to examine the distribution of radiation escaping from the scattering cloud in the case of the modified 
bremsstrahlung initial spectrum given by equation 16311 . We shall again assume that t csc = constant, and therefore we can 
employ equation l|82|l for the escape probability distribution P. In Figure 7, we plot the results for the escaping photon distri- 
bution N x (x, t) computed using equation 1)85^ . with f(x, y) set equal to the particular solution for the modified bremsstrahlung 
initial spectrum given by equation l|670 . Results are presented for two values of the low-energy cutoff x* and two values of 
the mean Compton parameter y. For small values of t — to, the escaping spectrum resembles the initial distribution. At inter- 
mediate times, the spectrum displays the hardening that is characteristic of Comptonized bremsstrahlung, and as t increases, 
the spectrum approaches the Wien distributution. The details of the spectral evolution depend on the values of x* and y, 
which may allow the extraction of these parameters from a sequence of spectra observed during a strong X-ray transient. For 
example, let us suppose that the transient begins with the impulsive injection into the corona of a pulse of soft bremsstrahlung 
radiation with an initial spectrum given by equation (1631 . As Figure 7 indicates, the shape of the spectrum for energies x < 1 
at the time of appearance of the Wien peak, combined with the overall amplitude of the spectrum at that time, is a sensitive 
indicator of the unique values of x* and y. In general, the low-energy spectrum increases with decreasing x* because it takes 
longer for the additional soft photons to be upscattered, and the overall amplitude of the spectrum at the time of appearance 
of the Wien peak increases with increasing y due to the competition between Comptonization and photon escape. Indeed, the 
simple observation of a Wien peak in the escaping spectrum would indicate that y>l, since otherwise the peak would occur 
after the normalization had already declined by several orders of magnitude due to photon escape, rendering the Wien bump 
unobservable. Increasing either x, or y results in more rapid equilibration to the Wien shape. The type of spectral evolution 
depicted in Figure 7 is expected to occur when the timescale for the impulsive injection of the bremsstrahlung radiation 
is shorter that the timescale for the radiation to escape from the cloud. In the next section we consider the situation with 
continual photon injection, resulting in the evolution of the spectrum towards a steady-state form. 

5.4 Generation of the steady-state Green's function 

The solution we have obtained for the Laplace transform L(x,xo,s) is closely related to the steady-state photon energy 
distribution resulting from the continual injection of monoenergetic radiation into a cloud of scattering electrons that is 
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Figure 7. Time-dependent escaping photon distribution N x (x, t) (eq. 1851 : solid lines) plotted in units of t^J for the case of the modified 
brcmsstrahlung initial spectrum fog. 1631 dotted lines). The elapsed time t — to in units of t CB c is indicated for each curve, and the values of 
the low-energy cutoff x* and the mean Compton parameter y are given by (a) x„ = IQ~ 3 , y = 0.2; (b) x t = 10 -3 , y = 1; (c) ie* = 10 — 2 , 
y = 0.2; and (d) x„ = 10~ 2 , y = 1. Included for comparison in each case is the Wien spectrum (filled circles) evaluated at the time tmax, 
where t m ax is the maximum time for the sequence of curves in the panel. 



leaking photons into the surrounding space. In this section, we explore the connection in detail, and establish that our time- 
dependent Green's function solution can be used to generate the steady-state Green's function corresponding to a balance 
between photon injection, thermal Comptonization, and photon escape. The steady-state solution is one of the most important 
and widely applied spectral models in X-ray data analysis, and it successfully accounts for the formation of the power-law 
spectra commonly observed in the quiescent emission from galactic and extragalactic sources. 

Let us suppose that the monoenergetic photon source is turned on at time t = 0. In this situation, the number spectrum 
of the photons remaining in the cloud at time t can be written in terms of the Green's function as the convolution 

N x (x,t) = ( N P(t)x 2 f G (x,x ,y)dt , (86) 
Jo 

where iVo expresses the rate of injection of fresh photons with dimensionless energy xo into the cloud per unit time, and 
y(t) = n e a T c(kT e /m e c 2 )(t — to) according to equation I1H . The number spectrum N x (x,t) is related to the total number of 
photons in the cloud N r (t) via N r (t) — N x (x,t) dx (see eq. 1761 ). Transforming the variable of integration in equation 1861 1 
from to to £ = (t — to)/t csc and substituting for P using equation I82H yields 

/•t/t esc 

N x (x,t)= te^Noe^x 2 f a (x,x ,y£)d£ , (87) 

Jo 

where y — t CBC n e a T c(kT e /m e c 2 ) (eq. (84j). The spectrum of the radiation escaping from the cloud is given by 

N x (x,t)= N e- £ x 2 f o (x,x ,y£)d4 . (88) 

Jo 

In the limit t — + oo, a steady-state balance is obtained between monoenergetic photon injection, thermal Comptonization, 
and photon escape. The resulting steady-state spectrum of the escaping radiation, A^ td (a;), is given by 



N?*(x)= Mm N*(x,t)= I N e- s x 2 f a (x,x ,y0dZ. (89) 

t — *oo 



By construction, this represents the Green's function solution for the time-independent problem. Transformation of the variable 
of integration from £ to y = y£ now yields 
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K td (x) = [°° e-*' s f a (x,x ,v)dy . (90) 

V Jo 

Comparing this expression with equation 1221 . we recognize that the integral on the right-hand side is simply the Laplace 
transformation of f G (x,xo,y) with the complex variable s replaced by 1/y. We can therefore immediately write 



1 

y \ y 



N* d (x)= 1 -^L[x,x ,~ I . (91) 



Substituting for L(x,xo, 1/y) using equation 125H yields the equivalent result 

K td (x) = y + Xo 2 e (x °~ x)/2 M 2 ,,(x min ) W^w) , (92) 
where a; ma x and a; m in are defined by equations 126H . and 
9 1^ 1/2 



This is essentially the same solution to the stationary Kompaneets equation derived by Sunyaev & Titarchuk (1980) in their 
Appendix B if we set their separation constant 7 = 1/y. The normalization of our solution for N^, td (x) can be confirmed by 
integrating equation 19211 to show that the total number of photons escaping from the cloud per unit time is given by 

N° td (x)dx = N , (94) 

as expected. Hence we have demonstrated that our solution for the time-dependent Green's function can be used to generate 
the steady-state Green's function. However, our solution also contains additional information describing the evolution of the 
spectrum towards the steady state, as discussed below. 



5.5 Evolution towards the steady-state spectrum 

We can use the results derived in § 5.4 to study the evolution of the photon distribution following the "turning on" of a 
continual source of monochromatic radiation at time t = 0. We assume that no photons are initially present in the cloud. 
This scenario describes, for example, an X-ray transient driven by the continual injection of iron line emission. Whether or 
not a steady state is established depends on whether the injection phase lasts longer than the characteristic timescale for the 
escape of radiation from the cloud. In Figure 8, we plot a sequence of curves describing the variation of the time-dependent 
escaping spectrum N x (x, t) computed using equation JHHJ, with f G (x,xo,y) evaluated using equation ifity . We also include 
the steady-state spectrum JV| td (a:) given by equation Ij92^ . Two different values of the initial photon energy xo and the mean 
Compton parameter y are considered. Note that as time increases, the spectrum broadens and the time-dependent curves 
approach the corresponding steady-state solutions as expected. The steady-state spectrum exhibits power-law wings around 
x — xo when y<l and xo < 1 (see Rybicki & Lightman 1979). For large values of x, the spectrum declines exponentially. 
The final equilibrium represents a balance between Comptonization, photon injection, and photon escape. The steady-state 
spectrum gives good agreement with the quiescent X-ray spectra observed from many active galaxies. However, with the 
availability of the new solution for the time-dependent Green's function, it is now possible to model the transient phase in 
the development of the steady-state spectrum using equation 18811 . This new theoretical framework improves our ability to 
interpret spectral/temporal data associated with flares that are driven by the injection of radiation on timescales shorter than 
that required for the establishment of a steady state. 



6 CONCLUSIONS 

In this paper, we have derived several new results of importance in the theory of time-dependent thermal Comptonization. 
The analysis is based on the transport equation (Qp, which describes the effects of energy diffusion (second-order Fermi 
energization), electron recoil, and spatial diffusion. The spatial and energetic components of the problem were separated by 
writing the Green's function for the occupation number n G as the product of the energy Green's function f G (x, xo,y) and the 
spatial distribution r/(r, fo, y) (eq. 1151 ). which is valid provided the scattering cloud has a steady, homogeneous structure. The 
Green's function n G represents the response of the system to the injection of a single photon with an arbitrary energy at an 
arbitrary time and location within the cloud. The particular solution for the occupation number distribution corresponding 
to a general source term can be obtained by integrating n G using equation J7J. 

The central result obtained in the subsequent analysis is the closed-form solution for the Kompaneets Green's function, 
f G (x, xo,y), given by equation 133H and plotted in Figure 2. This fundamental expression describes the evolution of an initially 
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Figure 8. Time-dependent escaping photon distribution N x (x,t) (eq. 1881 : solid lines) plotted in units of the injection rate TVo for the 
indicated values of the elapsed time t in units of t CBC . The scenario considered here corresponds to the continual injection of monochromatic 
radiation with energy xq, beginning at time t = 0. The values of xq and the mean Compton parameter y are given by (a) xo = 0.1, 
y = 0.2; (b) xq = 0.1, y = 1; (c) xq = 10, y = 0.2; and (d) xq = 10, y = 1. As time increases, the spectrum broadens and the 
time-dependent curves approach the steady-state solution N^ td (x) (eq. 1921 ; dashed lines). 



monoenergetic radiation spectrum in an infinite medium under the influence of repeated Compton scattering. As such, it serves 
as the "kernel" for thermal Comptonization, from which one can obtain the particular solution for any distribution of photon 
sources in time, space, and energy. The validity of this new analytical solution was confirmed in § 3 via comparison with the 
results obtained for the energy spectrum by integrating numerically the Kompaneets equation (see Fig. 2). The availability 
of our closed-form solution for the Green's function describing thermal Comptonization provides a valuable new tool for the 
analysis and interpretation of variable X-ray spectra. In particular, the Green's function represents an efficient alternative 
to numerical integration of the Kompaneets equation that avoids the complications and uncertainties associated with the 
imposition of boundary conditions in the energy space. 

The properties of the Green's function were explored in detail in § 4, where we confirmed that it displays proper equili- 
bration to the Wien spectrum at large times. We also reproduced Kompaneets' result concerning the variation of the mean 
photon energy x as a function of time (eq. 1461 : Fig. 3), once a typographical error made by Kompaneets is taken into ac- 
count. Furthermore, we obtained a new expression for the variation of the inverse-Compton temperature associated with a 
monoenergetic initial spectrum (Fig. 4), given by T IC /T e = (1/4) Zf//^, where the moments I§ and l2 are evaluated using 
equations 146H and 15211 . respectively. We also used our solution for the Green's function to derive a new, exact expression for 
the spectrum resulting from the time-dependent Comptonization of a bremsstrahlung initial spectrum, including a low-energy 
cutoff that approximates the effect of self-absorption. The new solution equilibrates to the Wien spectrum at large times (see 
eq. [S3 and Fig. 5). 

Although the Green's function f G (x, xo,y) obtained here represents the solution to the Kompaneets equation in an infinite 
medium, the formalism we have developed can also incorporate the effects of spatial transport in a homogeneous scattering 
cloud of any size and shape. In § 5, modifications associated with spatial transport in a cloud of finite size were examined 
using a simple escape-probability model. In this scenario, photons injected into the cloud with energy xo at some initial time 
t — to axe subject to thermal Comptonization as they diffuse out of the cloud. The analytical result for the escaping radiation 
spectrum N x (x,t) given by equation 18511 was plotted for monoenergetic and bremsstralung initial spectra in Figures 6 and 7, 
respectively. We also analyzed the case of the continual injection of monoenergetic photons into a "leaky" cloud, starting at 
time t = 0. A new time-dependent solution describing the gradual build-up of the escaping radiation spectrum was obtained 
(eq. |88j). and it was explicitly confirmed that in the limit t — > oo, the time-dependent spectrum approaches the steady-state 
solution given by Sunyaev & Titarchuk (1980), describing a balance between photon injection, thermal Comptonization, and 
photon escape (see eq. 1921 and Fig. 8). 
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Figure 9. Exact Green's function x 2 f G (x,xo,y) (eq. 1331 : solid lines) plotted as a function of the dimensionless photon energy x 
for the indicated values of the dimensionless time y. The initial photon energy is given by (a) xq = 0.1 and (b) xq = 1. Included for 
comparison is the "soft-photon Green's function," x 2 / (x, xq, y) (eq. 1981 : dashed lines), which ignores the recoil term in the Kompaneets 
equation 1161 . and therefore does not equilibrate to the Wien spectrum (1/2) x 2 e~ x (filled circles) as y — » oo. Conversely, f G (x,XQ,y) 
satisfies the complete Kompaneets equation, and therefore it approaches the Wien spectrum as y increases. 



6.1 Relation to previous solutions 

The new solution for the Green's function obtained in § 3.2 describes the fundamental physics of thermal Comptonization 
in a homogeneous plasma with steady properties. It therefore unifies and extends several of the previously known analytical 
solutions. For example, the new solution generalizes the analytical solution obtained by Zeldovich & Sunyaev (1969) and 
Payne (1980) for the Comptonization of soft photons with i< 1, which satisfies the equation 

9l = L±( x *0T\. (95) 
dy x 2 dx \ dxj 

This equation does not include the electron recoil term proportional to / that appears in the full Kompaneets equation 11611 . 
and therefore it only treats the stochastic energization of the photons. The exact solution to equation 195H . which we refer to 
as the "soft-photon Green's function," can be obtained using the Laplace inversion integral (eq. (29j ) by working in the limit 
i < 1, xq <C 1 (Titarchuk 2003, private communication). Along the inversion contour, Res = 7 > 0, and therefore it follows 
that Refj. > since fi = (s + 9/4) 1//2 according to equation |24| . The Whittaker functions Wb.^x) and M2,^(x) appearing in 
equation 1251 for the Laplace transform L(x,xo,s) consequently have the leading behaviors (see eqs. 1271 and 1281 ) 

M 2 „(x) ~ e~^x^ , W a ,„(z) = r{ l^ /2 f ~ X/2 ^ +1/2 > ( 96 ) 

for small x. Utilizing these asymptotic expressions for the Whittaker functions in equation 1251 . we find that the Laplace 
transform of the soft-photon Green's function is given by 

L solt (x,x ,s)= (xox y /2 (^Y , (97) 

where x max and Xmin are defined by equations 1261 . By performing the inverse Laplace transformation of Lsoft using equa- 
tion (1291 . we obtain for the soft-photon Green's function the solution 

(In x — Inxo) 2 



/soft (x,x ,y) = — exp 

2 jT^y 



which satisfies equation 195H . as well as the initial condition (see eq. 1171 ) 



(98) 



= Xq 2 S{x - x ) . (99) 



Note that due to the neglect of the recoil effect, / aoft fails to equilibrate to the Wien distribution as y — > 00. In Figure 9 we 
compare / aoft with the general Green's function / G (eq. 1331 ) for two values of the initial photon energy, xo = 0.1 and xo = 1. 
At early times (i.e., small values of y), the two solutions agree closely. However, as y increases, strong disagreement becomes 
evident. In particular, at high energies (large a;), the / so(t distribution continues to increase far beyond the saturation level 
indicated by the Wien spectrum. Conversely, our exact solution for / G , which includes the effect of electron recoil, equilibrates 
to the Wien form as required, and it therefore generalizes the Zeldovich & Sunyaev (1969) and Payne (1980) solution. 

In § 4.5 we considered the important problem of the Comptonization of a "modified" bremsstrahlung initial spectrum 
that includes a low-energy cutoff at energy x = x», which simulates the effect of self-absorption. Our new time-dependent 
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analytical solution for the resulting spectrum (eq. 1671 ) generalizes the analytical solution obtained by Chapline & Stevens 
(1973) and Becker & Begelman (1986), which sets the low-energy cutoff x» = 0, and therefore contains an infinite number of 
low-energy photons. While the x* = solution cannot equilibrate to the Wien spectrum, our new solution satisfies this physical 
requirement, and it therefore represents a distinct improvement over the previously known analytical result for bremsstrahlung 
injection. We demonstrated in § 4.5 that the low-energy cutoff x, and the mean value of the Compton y-parameter, y, can be 
estimated from observations of a sequence of spectra obtained during a bremsstrahlung-driven X-ray flare. Knowledge of these 
quantities yields constraints relating the electron temperature T e and optical thickness r of the corona to the gas temperature 
T and effective temperature T e s in the source region (see eqs. (621 and 1841 1. 

The new solution for the Green's function obtained here also generalizes the steady-state solution derived by Sunyaev 
& Titarchuk (1980) through the inclusion of the full time dependence of the spectral evolution process. However, it should 
be pointed out that the Sunyaev & Titarchuk formalism is also applicable to the case of an inhomogeneous scattering cloud, 
whereas the time-dependent solution obtained in the present paper was developed under the assumption of a homogeneous 
electron distribution. In future work we intend to investigate the possibility of generalizing the time-dependent Green's function 
in order to treat the inhomogeneous case as well. This would be a significant improvement since observations suggest that the 
electron distribution is probably inhomogeneous in a number of X-ray sources, as discussed in § 6.2. We briefly review the 
potential astrophysical relevance of the results presented in this paper below. 

6.2 Astrophysical relevance 

The closed-form solution for the Green's function describing thermal Comptonization (eq. 1331 1 provides a useful new tool 
for the analysis and interpretation of variable X-ray spectra. Analytical solutions offer far more physical insight than do 
numerical simulations, and they also provide the flexibility to efficiently explore a large region of parameter space. The 
new solution is relevant in a wide variety of astrophysical situations involving thermal Comptonization, such as studies of 
the reprocessing of the cosmological background radiation (Shimon & Rephaeli 2002; Colafrancesco et al. 1997), models of 
thermal upscattering in gamma-ray bursts (Ghisellini & Celotti 1999; Liang et al. 1999), studies of variability and lags in 
the X-ray spectra observed from AGNs and galactic black- hole candidates (Zdziarski et al. 2002; Lee et al. 2000), spectral 
reprocessing in X-ray bursts (Titarchuk 1988; Lapidus, Sunyaev, & Titarchuk 1987), and the Comptonization of line features 
in AGN spectra (Misra & Kembhavi 1998; Wang et al. 1999). The solution presented here is also perfectly suited for modeling 
the acceleration of relativistic electrons in a turbulent plasma. In this application, the terms proportional to df G /dx and f G 
inside the parentheses on the right-hand side of equation 11611 describe the second-order Fermi acceleration of the electrons 
and the effects of synchrotron/inverse-Compton losses, respectively (Schlickeiser, Sievers, & Thiemann 1987; Schlickeiser 1984; 
Lacombe 1977). The availability of our exact solution for the Green's function therefore provides a unique opportunity to 
model analytically the time-dependent development of the electron energy spectrum during solar flares and other high-energy 
astrophysical transients (e.g., Park, Petrosian, & Schwartz 1997; Miller & Ramaty 1989). 

The best prospects for extragalactic, time-resolved X-ray spectroscopy of accretion flows around black holes are presented 
by nearby AGNs such as NGC 4051, NGC 4151, NGC 6814, and MCG-6-30-15, which are among the brightest and most 
strongly variable X-ray sources known. The rapid variability of these and other Seyfert galaxies represents one of the most 
interesting puzzles in modern X-ray astronomy. A large body of observational evidence suggests that as the 2-10 keV flux 
increases in these sources, the X-ray spectrum becomes softer (Lee et al. 2000; Merloni & Fabian 2001; Nandra 2001). For 
example, observations of the dwarf Seyfert nucleus of NGC 4395 reported by Iwasawa et al. (2000) reveal rapid variability, 
including flares during which the flux changes by a factor of 3-4 in a single half-day ASCA observation. Doubling times as 
short as 100 s have been observed. The highly variable Seyfert galaxy NGC 6814 has displayed factor of 3 changes in flux 
in 8 hours (Konig et al. 1997). The characteristic softening observed during the Seyfert flares may reflect inverse- Compton 
cooling of the corona due to the injection of copious soft photons (possibly due to bremsstrahlung) from the inner region 
of the accretion disk. The spectral evolution observed during these transients may provide us with a direct opportunity to 
study time-dependent Comptonization. However, it is also possible that the variable X-ray spectra may be more correctly 
interpreted using a sequence of steady-state models, if the timescale for the injection of the photons exceeds the characteristic 
time for the photons to diffuse out of the cloud (Sunyaev & Titarchuk 1980). 

The signature of time-dependent Comptonization is perhaps most unambiguously observed in the complex cross-spectra 
produced by combining the Fourier transformed lightcurves at two different X-ray energies (van der Klis et al. 1987). The 
resulting structure of the X-ray time lags contains detailed information on the "Compton reverberations" taking place in the 
inner region of the disk/corona (Payne 1980; Reynolds et al. 1999) that cannot be obtained using standard spectral analysis 
(Kazanas, Hua, & Titarchuk 1997). Observations of hard X-ray time lags in AGNs and galactic black-hole candidates are 
consistent with the predictions of thermal Comptonization models (Zdziarski et al. 2002; Lee et al. 2000). The time lag is 
a general consequence of the upscattering of soft photons by hot electrons (Payne 1980). It is worth emphasizing that the 
upscattering of soft photons can take place even in situations where the electron temperature is decreasing as a function of 
time due to inverse-Compton cooling, as discussed in § 6.3. Hence the spectral softening, combined with the hard lags, strongly 
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imply that thermal Comptonization is playing a significant role in these sources. The hard lags are a particularly important 
diagnostic, since they may provide an indelible "imprint" of the time-dependent process that can be detected even when there 
is inadequate spectral-temporal resolution to accurately measure the sequence of spectral curves. 

The time-dependent nature of the analytical solution for the Green's function obtained here provides a convenient basis 
for the theoretical computation of the frequency-dependent lags. By exploiting the close relationship between the Fourier and 
Laplace transforms, one can use the exact expression for the Laplace transformation of the energy distribution (eq. |25| ~) along 
with the Laplace transformation of the escape probability distribution to obtain the Fourier transformation of the escaping 
spectrum using the convolution theorem. This would yield an analytical means for computing the complex cross-spectrum 
and evaluating the X-ray time lags (van der Klis et al. 1987), which may represent an efficient alternative to the Monte-Carlo 
based simulations that are usually employed to compute the lags. The formalism developed here focuses on a homogeneous 
scattering cloud, and there are some indications that this may not be an appropriate model for certain galactic black-hole 
candidates, in which the density in the corona may vary as 1/r (Bottcher & Liang 1998; Kazanas, Hua, & Titarchuk 1997). 
However, this conclusion also depends on the shape of the cloud, and it would therefore be interesting to extract the lags 
using the analytical solution and compare them with the available X-ray data. Moreover, it may be possible to extend the 
general analytical approach taken here in order to treat inhomogeneous clouds. We plan to pursue this possibility in future 
work. 

In addition to the continuum emission discussed above, Seyfert galaxies also appear to exhibit Fe Ka line emission. This 
component of the spectrum is quite broad, and in some cases highly variable (Wang, Zhou, & Wang 1999; Wang et al. 1999; 
Vaughan & Edelson 2001). It is currently not clear whether the width of the line features reflects Doppler and gravitational 
broadening very close to the black hole (Reynolds & Wilms 2000), or perhaps broadening due to thermal Comptonization in 
a surrounding corona (Misra & Kembhavi 1998). Our analytical solution for the Green's function may prove to be helpful 
in resolving this issue since it provides an exact description of the time-dependent reprocessing of the line emission in the 
corona. Sample calculations illustrating the broadening of the line were presented in Figure 6, where the effects of variations 
in the initial energy eo relative to the electron thermal energy kT e were explored, along with variations in the value of the 
mean Compton parameter y. The analytical expression for the time-dependent Green's function developed here is valid for 
arbitrary values of the initial photon energy eo, and therefore it can be used to study the upscattering of an initial photon 
distribution with eo <§C kT e , or the downscattering of initial photons with eo 3> kT e . We reiterate that the previous analytical 
solution for the time-dependent Green's function derived by Payne (1980) and Zeldovich & Sunyaev (1969) is valid only for 
the case of very soft photons scattered by hot electrons, due to the neglect of the recoil term in the Kompaneets equation. 

6.3 Thermal and dynamical effects 

It has long been recognized that the injection of large quantities of soft radiation into a hot corona will almost certainly result in 
substantial inverse-Compton cooling of the electrons (Payne 1980; Lightman, Giacconi, & Tananbaum 1978; Guilbert, Fabian, 
& Ross 1982; Becker & Begelman 1986). This idea is further supported by the fact that the spectral softening observed in 
Seyfert X-ray transients is inconsistent with the hardening evident in the sequence of Comptonized bremsstrahlung spectra 
plotted in Figure 7, based on the assumption of a constant-temperature corona. We therefore conclude that the variation of 
the electron temperature in the corona must be included in models for the spectral evolution occurring during the flare. If 
the transient is very intense, the gas is expected to be radiation-dominated, and the electron temperature should therefore 
track the inverse-Compton temperature of the radiation (Becker & Begelman 1986; Becker 1999). This issue has been explored 
by Bottcher (2001), who performed Monte-Carlo simulations of disk/corona models with temperatures that responded self- 
consistently to the inverse-Compton cooling. They find that the inclusion of thermal effects improved the agreement between 
the predicted and observed X-ray time lags. If the cooling of the electrons occurs on timescales that exceed the characteristic 
time for the photons to diffuse out of the cloud, then the spectrum may be adequately represented using a sequence of 
constant-temperature (though not necessarily time-independent) models. On the other hand, if the electron temperature 
varies on shorter timescales, then the entire problem of thermal Comptonization must be reconsidered from first principles, 
including a self-consistent treatment of the temperature. It is currently unclear whether the temperature varies this rapidly 
during typical Seyfert transients (see discussions in Payne 1980; Petrucci et al. 2001; Nandra 2001; Malzac & Jourdain 2000). 
We plan to explore this question in future work. 

Another issue of potential importance for models of the spectral evolution during X-ray transients is the possible dynamical 
response of the corona, driven by the heating and cooling of the gas. Bulk motions may be important if the gas responds 
hydrodynamically to the energy input, or if it is already moving, as in an accretion flow. Expansion of the corona due to 
heating may deplete some of the energy from the photons, and, conversely, contraction will lead to photon energization via 
the first-order Fermi process ("bulk Comptonization"). Although we have assumed here that the electron temperature T e and 
number density n e remain constant during the transient, in general these two quantities should be considered functions of 
time and position, since both will respond to the hydrodynamics and thermodynamics of the gas. Some progress has been 
made in understanding the effect of the bulk flow on the radiation spectrum in idealized situations. For example, Colpi (1988) 
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has analyzed the effect of thermal Comptonization on photons injected at an arbitrary radius into a freely-falling, spherical 
accretion flow using a steady-state model. She finds that the convergence of the flow can significantly alter the shape of 
the escaping high-energy spectrum due to a combination of photon trapping and Fermi energization. Similarly, Laurent & 
Titarchuk (2001) find that bulk Comptonization in quasi-spherical inflows can have a direct effect on the steady-state power- 
law continuum emission produced. Following a similar approach, it may be possible to incorporate the effects of bulk motion 
within the framework of a time-dependent model such as the one considered here. 

In conclusion, we believe that the theoretical framework developed in this paper will facilitate a more complete utilization 
of the high-quality temporal and spectral data provided by current and future orbital X-ray observatories. The new solution 
for the Green's function generalizes and extends the previously known analytical solutions for thermal Comptonization, 
and further extensions to cases involving variable coronal temperatures and density structures may also be possible. When 
combined with time-resolved spectra obtained during bright X-ray flares in AGNs, these analytical solutions may allow the 
determination of source parameters that previously were difficult to interpret due to the necessity of running lengthy computer 
simulations. 

The author is grateful to Demos Kazanas for several stimulating conversations, and also to the referee, Lev Titarchuk, 
for providing many useful comments which led to significant improvements in the manuscript. 
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APPENDIX A: MATHEMATICAL DETAILS 
Al Formal solution for the Laplace transform 

The Laplace transform 

L(x,x ,s) = / e~ sy f G (x,x ,y)dy (Al) 
Jo 

of the Green's function distribution f G (x,xo,y) satisfies the inhomogeneous ordinary differential equation 



( L+ dL\ 
\ dx 



s L — —x 8(x — xo) ■ (A2) 



l_d_ 

x 2 dx 

The homogeneous equation obtained when x 7^ xq has the solutions 

!A x~ 2 e~ x/2 M 2:ll (x) , x < x , 
(A3) 
B x- 2 e~ x/2 W 2 ,„(x) , x > x , 

where M 2tll {x) and W 2 ^(x) denote Whittaker's functions, and 

1/2 



/' 



/9 \ 1 ' 2 

(*)=(-+<) • (A4) 

Asymptotic analysis indicates that these are the only functions that yield convergent results for the photon number and 
energy densities. The constants A and B appearing in equation (IA3I are determined by imposing continuity and derivative 
jump conditions on the transform L. The solution for L must be continuous at x = xo in order to avoid an infinite flux of 
photons in the energy space. The continuity condition requires 

AM 2 , ll (x ) = BW 2 , li (x ) . (A5) 
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Furthermore, the derivative dL/dx must display a jump at x — xo given by 



dL 
lim — — 

e->o dx 



Xq+E 



dL 
dx 



which is obtained by integrating equation lAH in a narrow range around x — xq. 

Utilization of the continuity and derivative jump conditions yields for A and B the solutions 

■e xo/2 Xo 2 W2,n(x ) D -e xo/2 x- 2 M2, l ,{x ) 

ID — 



A 



w(xo) w(xo) 
where the Wronskian of the two solutions is defined by 

w(x) = M 2 , »{x) ^- W 2 , - W 2 ,„(x) ^- M 2 ,„(x) . 



(A6) 



(A7) 



(A8) 



The Wronskian can be evaluated using equations (13.1.22), (13.1.32), and (13.1.33) of Abramowitz & Stegun (1970), which 
yields 

r(l + 2/i) 



w{x) - TQ* -a/2)' 

Combining results, we find that the formal solution for the Laplace transform is given by 

L(x,x ,s) = r( ^- 3 ^ x~ 2 x- 2 e^ 2 

K N lM 2 ^{x )W 2 ^(x) , x>x , 

or, equivalently, 



L{x, xo, s) 



rfr-3/2) a a (.„■ 
r(i + 2 M ) Xo 



r)/2 



M 2 , M (a; m in) W 2 , M (a; max J 



where 

i min = min(x,x ) 



max(i, xo) 



(A9) 



(A10) 



(AH) 



(A12) 



The solution for the Laplace transform L(x,xo,s) is closely related to the Green's function for the steady-state problem 
developed by Sunyaev & Titarchuk (1980), as discussed in § 5.4. Simple poles are located where the function F([i — 3/2) 
diverges, which occurs when the quantity fi — 3/2 is zero or a negative integer. The evaluation of the residues at the simple 
poles is carried out in the next section, and the inverse Laplace transformation is discussed in § 3.2. 



A2 Evaluation of the residues 

We found in § 3 that the Green's function energy distribution can be expressed as 



f G (x,x ,y) 



2ni 



e sy L(x, xo, s) ds - 



1 f Q 2 
-— / e ay L(x,xo,s)ds+ > Res(s„) , 
2m J p ^ 



(A13) 



in the limit n — ► oo, r 2 — > (see Fig. 1), where the transform L(x, xo, s) is given by equation IA10I and Res(s n ) is the residue 
associated with the pole at s = s n . In our problem, simple poles are located at si = and s 2 — —2, and the corresponding 
residues are computed using the formulas (Butkov 1968) 



Res(si) = lim sL(x,xo,s)e s 



Res(s2) = lim (s + 2) L(x, xo, s) e s 

a^- 2 



(A14) 



Because the poles correspond to the singularities of the function F(/i — 3/2), evaluation of the residues requires utilization of 
the asymptotic results 

lim sT(u-3/2) =3 , lim (s + 2) T(u - 3/2) = -1 , (A15) 

which follow from the definition of fj,(s) given by equation ljA4fl . We can show that the functions M 2 ^(x) and W 2yM ,(x) 
reduce to combinations of Laguerre polynomials and exponentials at the points s = and s = —2. In particular, using 
equations (13.1.2), (13.1.3), (13.1.32), and (13.1.33) from Abramowitz & Stegun (1970), we obtain for s = the results 



M 2 ^(x) 
Likewise, for s 



-x/2 2 

e x 



2 we find that 

1/2 x (l-|) , W 2 ,,(x) 



-x/2 2 

e ' x 



- x/2 x{2-x) 



(A16) 



(A17) 
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Combining equation 1A101 with equations JA14I - II A 1 71) . we establish that the residues are given by 



Res(si) = - e 



Res(s2) = - e 



-x-2y 



(2 — x) (2 — Xq) X X 



(A18) 



A3 Integration along the branch cut 

The remaining integrations in equation IIA13fl along the upper and lower segments NO and PQ, respectively, need to be 
handled carefully because of the presence of the branch cut extending from s — —9/4 to s — — oo (see Fig. 1) in the definition 
of /x(s) (eq. IA4I 1. We proceed by transforming to a new variable of integration, it, defined by 



9 

"4" S 



(A19) 



with u = at the branch point. This is useful because the integration bounds in equation 1A13II now transform to the 
convenient values (0, oo). In terms of it, equation l|A40 for u(s) can be rewritten along the upper and lower segments as 



At (it) = i u , segment NO 
/x(it) = —iu , segment PQ 



(A20) 
(A21) 



The sign change between the two expressions for /Lt(u) along the two segments reflect the two different branches of the complex 
square root function. Using equation 1A111 to substitute for L{x,Xq,s) in equation 1AI3II and transforming the variable of 
integration from s to it, with ds = —2udu, we obtain 



f G (x,x ,y) 



r(-3/2- w) 

r(i - 2iu) 



M 2 , 



i(£min) W 2 , -iu (a; max I 



r(-3/2 + m) 



M 2 , in (Xuun) W 2 , iu (X n 



i du + S Res(sn) 



r(l + 2iu) 

where the residues are given by equations lAlSfl . From the definition of the W2,^{x) function, 



W 2 ,„(x) 



r (- 2 M) 1, ^^ , r ( 2 A0 



r(-/x-3/2) 

it follows that 



M 2:fl (x) + 



F(n - 3/2) 



M; 



2, —a 



(x) 



W2, iu (x) =W 2 ,-iu(x) , 

and therefore we can rewrite our expression for the Green's function as 

"r(-3/2 - iu 



(A22) 



(A23) 



(A24) 



f G (x,x ,y) 



W 2 , iu(x max ) 



F(I - 2iu) 



Ma, -™(a;miii) 



r(-3/2 + iu) 



M2,iu(Xjniu) 



i du + S Res(s n ) 



r(l + 2iu) 

Based on equation 1A23L we can show that the factor in square brackets in equation (IA25I is given by 

r(-3/2-m)r(-3/2 + in), 



r(-3/2-m) , , r(-3/2 + i«) , 
— : - -M 2 ,_TO(o; m in) . / M 2yIU (x lnin ) 



rri - 2iu) 



T(l + 2iu) 



T(2m) r(l - 2m) 



-W 2 , iu{Xmiii) 



The products of gamma functions on the right-hand side of this expression can be simplified using the identities 



T(2iu) r(l - 2iu) = 



2i sinh(7ru) cosh(7ru) 



I6tt 



r(-3/2 - iu) T(-3/2 + iu) = — r— , , , , , ., . 

cosh(7ra) (I + 4u 2 ) (9 + 4m 2 ) 

Combining equations JXT5|, JX25J, llA26l . SK571 . and llA28l . we obtain after some algebra 



(A25) 

(A26) 

(A27) 
(A28) 



f G (x,x ,y) 



32 



-9y/4 



2 (x -x)/2 



u sinh(7rit) 



(I + 4u 2 )(9 + 4ii 2 ) 



w , m„ , u , e~ x , e""" 22 ' (2 - a;) (2 - xo) 
x Wj,i«(xo) W2,i«(a:)du + — + 



(A29) 



2 so a; 

This is our final result for the Green's function, which appears as equation l|33fl in the main text. Series expansions for the 
Whittaker functions are provided in § A. 4. A self-contained FORTRAN program that evaluates the Whittaker functions and 
performs the integration in equation <A29t to yield the Green's function is available from the author upon request. We point 
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out that a related integral involving a single Whittaker function was worked out by Kompaneets (1957) in collaboration with 
I. M. Gelfand in order to describe the variation of the mean photon energy due to thermal Comptonization, as discussed in 
§ 4.2 (see eq. \M ). 



A4 Series expansions for the Whittaker functions 

In order to carry out the integration in the analytical solution for the Green's function given by equations 1331 and 1A29L we 
need to be able to evaluate the Whittaker function W2,m(x) for both large and small values of the argument x. Furthermore, 
the solution we have obtained for the optically thin bremsstrahlung initial spectrum (eq. |67p requires the evaluation of the 
functions Wo,iu(x*), W-i, iu(x*), and W-i,m{x*) for small values of the low-energy cutoff a;* (see eq. E]). We 

shall develop the necessary expansions below. 



A4-1 Expansion for Small x 

Based on the definition of W Klll (x) given by equation (13.1.34) of Abramowitz & Stegun (1970), we can write 

w kM *) s r(1/ r 2 ( : 2 ;^ m) "«.«.(*) + M — <*> • (A30) 

The two terms on the right-hand side are complex conjugates of each other, and therefore W K ,iu(x) is purely real. Hence we 
obtain 

W K , lu {x) = 2Re/ .T ( ~ 2m) . r M K , iu (x)\ . (A31) 



r(i/2 -K-iu) 

For small values of x, the function M Ktll (x) can be evaluated using the series 

M K: iu (x) = e"* /2 x^ Y (1/ ''" + m) " ^ , (A32) 

t—' (1 + 2lU) n TV. 
n— 

where (a) n denotes the Pochhammer symbol, defined by (Abramowitz & Stegun 1970) 

(o)„ = i ' . . . ' (A33) 

I a(a + 1) ■ ■ ■ (a + n — 1) , n > 1 . 

Combining expressions, we obtain finally 

W ^)=2e^^±^Bsl j^K , • (A34) 

^ n! 1 r(l/2 — k — iu) (1 + 2iu) n \ 

n— 

This expansion can be used to evaluate each of the Whittaker functions Wx,i u {x*), Wo,i U (x,), W-i t i U (x,), and W-2. iu(x*) 
appearing in equation 1671 . 

An application of special interest is the function W2,iu(x) appearing in equation il'-i-'il for the Green's function. We obtain 
in this case 

W 2 iu (x) = 2 S<> f ^ Re ( ff^ ^ 2 + ) . (A35) 
' w ^ n! lr -3/2-w (l + 2m) n j v y 

In practice, we use this expansion to evaluate W2, iu(x) for x < Ximd(u), where 

( \ - J 20 ' «< 10/3 , 

Xmid(M) = (lO + 3u, u>10/3 . (A36) 

We emphasize that the functions W K , iu(x) and W2, iuix) are purely real, despite the appearance of the imaginary number i 
in equations IA34II and IA35I I. The complex gamma function F(z) can be evaluated using equation (6.1.48) from Abramowitz 
& Stegun (1970). 



A4-2 Expansion for Large x 

The Whittaker functions Wi,i u (x*), Wo,iu(x,), W—i t iu(x*), and W-2,iu(x,) appearing in equation 1671 for the bremsstrahlung 
particular solution need only be evaluated for small values of the low-energy cutoff x*, and therefore equation l|A340 is sufficient 
for this case. However, the function W2, iu(x) appearing in equation l|33|l for the Green's function must be evaluated for general, 
positive values of x. We can use equation l|A350 to compute Wz,iu(x) when x < x m id(u). For larger values of x, a convergent 
series for W 2 ,iu{x) can be obtained by combining equations (13.1.33) and (13.5.2) of Abramowitz & Stegun (1970), which 
yields 
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oo 



(-3/2 + w)„(-3/2 



E 



— lU) n 



(A37) 



The Pochhammer symbols in this equation are complex conjugates of each other, and therefore our result for W2, j«(a;) can 
be rewritten as the completely real expression 



We use this formula to evaluate W2, iu{x) when x > x m id(u), where imid (it) is defined by equation l|A360 . In performing the 
numerical integration in equation I33L it is sufficient to treat the domain < u < 20 because larger values of u make a 
negligible contribution to the integral. 

This paper has been typeset from a TprjX/ I^Tf^X file prepared by the author. 
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